## Warning: Removed 18659 rows containing missing values (geom_path).
Look at raw data
## Warning: Removed 35 rows containing missing values (geom_path).
## Warning: Removed 16 rows containing missing values (geom_path).
How much data is missing?
** A lot of data missing in April let’s throw that out
## Warning: Removed 65 rows containing non-finite values (stat_boxplot).
Look at seasonal monthly median
## Warning: Factor `month` contains implicit NA, consider using
## `forcats::fct_explicit_na`
Seasonal monthly mean
## Warning: Factor `month` contains implicit NA, consider using
## `forcats::fct_explicit_na`
## Warning: Removed 3 rows containing missing values (position_stack).
## Warning: position_dodge requires non-overlapping x intervals
## Warning: Removed 18027 rows containing missing values (geom_col).
## Warning: Removed 2192 rows containing missing values (geom_path).
## Warning: Removed 14206 rows containing missing values (position_stack).
## Scale for 'y' is already present. Adding another scale for 'y', which
## will replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which
## will replace the existing scale.
## Warning: Removed 3 rows containing missing values (geom_point).
# SNOTEL and PRSIM
bs_snotel_prism <- bs_sc_snotel_prism_snodas %>%
mutate(month = month(date, label = TRUE, abbr = TRUE)) %>%
filter(month %in% c("Nov", "Dec", "Jan", "Feb", "Mar")) %>%
select(date, month, SNOTEL, PRISM) %>%
drop_na() %>%
mutate(Date = date) %>%
addWaterYear()
graph_correlation(bs_snotel_prism, x_data = bs_snotel_prism$PRISM,
y_data = bs_snotel_prism$SNOTEL,
title = "Beaver Spring SNOTEL Daily SWE and PRISM ",
xlabel = "PRISM Accumulated Snowfall (mm)",
ylabel = "SNOTEL SWE (mm)") +
scale_y_continuous(expand = c(0,0), limits = c(0,400))+
scale_x_continuous(expand = c(0,0), limits = c(0,400))+
facet_wrap(~waterYear)
## Scale for 'y' is already present. Adding another scale for 'y', which
## will replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which
## will replace the existing scale.
## Warning: Removed 45 rows containing missing values (geom_point).
# SNOTEL and SNODAS
bs_snotel_snodas <- bs_sc_snotel_prism_snodas %>%
mutate(month = month(date, label = TRUE, abbr = TRUE)) %>%
filter(month %in% c("Nov", "Dec", "Jan", "Feb", "Mar")) %>%
select(date, month, SNOTEL, SNODAS) %>%
drop_na() %>%
mutate(Date = date) %>%
addWaterYear()
graph_correlation(bs_snotel_snodas, x_data = bs_snotel_snodas$SNODAS,
y_data = bs_snotel_snodas$SNOTEL,
title = "Beaver Spring SNOTEL Daily SWE and SNODAS ",
xlabel = "SNODAS SWE (mm)",
ylabel = "SNOTEL SWE (mm)") +
scale_y_continuous(expand = c(0,0), limits = c(0,500))+
scale_x_continuous(expand = c(0,0), limits = c(0,500))+
facet_wrap(~waterYear)
## Scale for 'y' is already present. Adding another scale for 'y', which
## will replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which
## will replace the existing scale.
## Warning: Removed 20 rows containing missing values (geom_point).
## Error in as.character(x): cannot coerce type 'closure' to vector of type 'character'
## Warning: Removed 14155 rows containing missing values (position_stack).
## Warning: position_dodge requires non-overlapping x intervals
## Warning: Removed 17976 rows containing missing values (geom_col).
## Warning: Removed 2051 rows containing missing values (geom_path).
## Scale for 'x' is already present. Adding another scale for 'x', which
## will replace the existing scale.
## Scale for 'y' is already present. Adding another scale for 'y', which
## will replace the existing scale.
## Warning: Removed 2 rows containing missing values (geom_point).
## Scale for 'colour' is already present. Adding another scale for
## 'colour', which will replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which
## will replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which
## will replace the existing scale.
## Scale for 'y' is already present. Adding another scale for 'y', which
## will replace the existing scale.
## Scale for 'y' is already present. Adding another scale for 'y', which
## will replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which
## will replace the existing scale.
## Scale for 'y' is already present. Adding another scale for 'y', which
## will replace the existing scale.
## Scale for 'y' is already present. Adding another scale for 'y', which
## will replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which
## will replace the existing scale.
## Scale for 'y' is already present. Adding another scale for 'y', which
## will replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which
## will replace the existing scale.
## Warning: Removed 2 rows containing missing values (geom_point).
## Scale for 'y' is already present. Adding another scale for 'y', which
## will replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which
## will replace the existing scale.
## Scale for 'y' is already present. Adding another scale for 'y', which
## will replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which
## will replace the existing scale.
## Scale for 'y' is already present. Adding another scale for 'y', which
## will replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which
## will replace the existing scale.
## Warning: position_dodge requires non-overlapping x intervals
## Warning: Removed 5561 rows containing missing values (geom_col).
## Scale for 'y' is already present. Adding another scale for 'y', which
## will replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which
## will replace the existing scale.
## Warning: Removed 17 rows containing missing values (geom_point).
## Warning: position_dodge requires non-overlapping x intervals
## Warning: Removed 5558 rows containing missing values (geom_col).